Hall response of interacting bosonic atoms in strong gauge fields: from condensed to fractional 

quantum Hall states. 
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Interacting bosonic atoms under strong gauge fields undergo a series of phase transitions that take the cloud 
from a simple Bose-Einstein condensate all the way to a family of fractional-quantum-Hall-type states [M. 
Popp, B. Paredes, and J. I. Cirac, Phys. Rev. A 70, 053612 (2004)]. In this work we demonstrate that the 
Hall response of the atoms can be used to locate the phase transitions and characterize the ground state of the 
many-body state. Moreover, the same response function reveals within some regions of the parameter space, 
the structure of the spectrum and the allowed transitions to excited states. We verify numerically these ideas 
using exact diagonalization for a small number of atoms, and provide an experimental protocol to implement 
the gauge fields and probe the linear response using a periodically driven optical lattice. Finally, we discuss our 
theoretical results in relation to recent experiments with condensates in artificial magnetic fields [ L. J. LeBlanc, 
K. Jimenez-Garcia, R. A. Williams, M. C. Beeler, A. R. Perry, W. D. Phillips, and I. B. Spielman, Proc. Natl. 
Acad. Sci. USA 109, 1081 1 (2012)] and we analyze the role played by vortex states in the Hall response. 

PACS numbers: 03.75.Hh, ()3.75.Kk, 67.40. Vs 



I. INTRODUCTION 

Within the last few years it has been demonstrated that cold 
gases of neutral atoms play an important role as quantum emu- 
lators, as these systems provide controllable devices available 
to simulate other systems of interest or theoretical models. 
The great flexibility of cold gases for which different types 
of atoms or interactions as well as their environment can be 
selected nearly at will, offer many possibilities 1 1-3]. Among 
them is the simulation of the quantum Hall effect. 

So far two main alternatives to engineer sufficiently strong 
artificial gauge fields to generate fractional-quantum-Hall 
(FQH) states have been proposed. One is by atom-laser cou- 
pling and the other using rotating potential traps. When a neu- 
tral atom moves in a properly designed laser field, its center- 
of-mass motion may mimic the dynamics of a charged particle 
in a magnetic field. When the atom follows adiabatically one 
of its dressed states (i.e., local eigenstates of the atom-light 
coupling), artificial magnetism emerges, due to the accumu- 
lation of the Berry phase |jll 01- An alternative procedure to 
generate gauge fields is by the rotation of the trap potential 
that confines the system [5, 6], where the rotation frequency 
Q, plays the role of the magnetic field. In the regime of high 
magnetic flux, or rapid rotation in the case of a bosonic cloud, 
theory predicts the appearance of strongly correlated phases. 
These phases can be viewed as the bosonic version of FQH 
states. IVIoreover, many-particle systems provide a rich vari- 
ety of different phases of quantum matter. Different proposals 
of classification can be found in literature 

Nowadays it is experimentally feasible to create artificial 
fields for ultracold atoms trapped in optical lattices (see ifToll . 
Chap. 4 where a rather exhaustive analysis is made of the 
state of the art). Proposals for Abelian and even non-Abelian 
gauge potentials exist, which in some cases are candidates to 
exhibit Hall effect under special setups lfTTl - [T6ll . and further- 



more, some recent articles refer to valuable schemes to mea- 
sure signatures of the Hall effect [17-20|. A new generation 
of proposals to generate FQH states is given by Refs ll2lll22ll . 
However, experimental difficulties have prevented the obser- 
vation of strongly correlated states. On the one hand, the use 
of laser beams comes with some drawbacks, such as heating 
of the atoms due to residual spontaneous emission |1], and 
on the other hand, in the case of rotating trap, if the number 
of atoms is large (as is usually the case in experiments), the 
TV-dependent critical rotation frequency needed to enter the 
strongly correlated regime is so close to the trap frequency that 
the system becomes unstable. Nonetheless, vortex states have 
been experimentally obtained using both techniques [23, 24]. 

In practice, the measurement of interesting observables in 
an experiment with few atoms and high density, the regime 
in which interactions are most relevant, is very challenging 
because we may not have access to the spatial density pro- 
file. With this in mind, we present an experimental proposal 
to simulate the few-particle system in each site of an opti- 
cal lattice, wherewe choose the rotating trap-potential alterna- 
tive. Several interesting features are stressed: First, it is an 
ideal playground to test the properties of a Hall system with 
strong interactions. Second, the suggested experiment allows 
the study of a larger number of atoms than what is computa- 
tionally feasible with our techniques. Third, as the window 
of observability (the range of Q, values) of strongly correlated 
states increases as N decreases, the critical value of Q, from 
which FQH states of few atoms are observable is far from the 
instability region. Or in other words, taking advantage of the 
possibility to have few atoms per site {N < 10) we solve the 
experimental difficulty that has so far prevented the observa- 
tion of FQH states for large systems. And finally, we introduce 
a method to measure the linear response using time-of-flight 
images, which do not require spatial resolution of the atomic 
profile in each lattice site. 



2 



Measuring characteristic properties of the highly correlated 
state is a very relevant problem. We will show below that 
the linear response function related to the Hall effect provides 
valuable information about the eigenstates and phase transi- 
tions in the strongly correlated atomic system. This requires a 
measurement of the population of the "scissor mode" ({xy)), 
which is Hall excited by the "breathing-mode" perturbation 
Hpert cos{ujt); see |25] for the first experimental obser- 
vation of the Hall effect with atoms in the mean-field regime. 
We expect that experimentally the same informationwill be 
available in the Hall response of the system, at least in the 
weak perturbation limit. We use the rotating frequency as 
the driven parameter and find that as the rotation increases, 
new phases emerge which are directly related to angular mo- 
mentum transitions |26-'29'l. To characterize the nature of the 
many-body ground states within a phase, we analyze the role 
played by its excitation spectrum in the Hall dynamical re- 
sponse. We find that the Hall response increases at phase tran- 
sition points, and that it is modified in the presence of vortex 
states. 

This paper is organized as follows: In Sec. II we present the 
model of the unperturbed system in two parts. In Sec. IIA we 
show the equations used to perform exact diagonalization of a 
rotating system to obtain the ground state (GS) and the excita- 
tions. In Sec. IIB we propose an optical lattice setup with the 
implementation of rotating traps in each independent site and 
infer the expression of the trap frequency used in Sec. IIA as 
a function of the experimental parameters. In Sec. Ill we dis- 
play the expressions used within the linear response theory to 
obtain the Hall response function at an individual site and how 
it can be implemented in the lattice. In Sec. IV we analyze 
our results and give an interpretation. Finally in Sec. V we 
present our conclusions. In the Appendix, we give a detailed 
explanation of our optical lattice proposal. 



n. MODEL OF THE UNPERTURBED HALL SYSTEM 

A. Analytical background: A single rotating trap 

Inside each site of the optical lattice, we assume a twodi- 
mensional system of N bosonic atoms of mass M. The cloud 
is trapped in a rotating parabolic potential of frequency uj± 
and rotation fl along the z axis. In the rotating frame of refer- 
ence, the Hamiltonian reads ll30ll 



Ho — Hsp+Hint, (1) 

the single particle (sp) part given by, 
with 



where the particular selection of the symmetric gauge has 
been made in the definition of A, being B* = 2M^ a con- 
stant artificial magnetic field directed downward along the z 
direction and r — {x, y). From now on we consider M = 1/2 

and h = 1 and choose — \J aiIj^ ~ ■\/2/ajj_ , fujj±/2 
and as units of length, energy, and frequency, respec- 

tively. With our unit of length, wj^ = 2. W fixes a term that 
breaks the isotropy of the trapping potential and is given by 

W{x,y)^2TMu:l{x^ -y^) (4) 

where the dimensionless parameter r measures the strength 
of the anisotropy. With this term, the part of the trapping po- 
tential which is independent of B* can be rewritten as 

Vtrap{x,y) = ^(w^a;^ + Wyy^) (5) 

where = + 4t) and uj^ = — 4r), being 

r < 1/4. For the sake of stability we require B*) < 
u}j_y/l ~ 4t = 2y/l - At. 

We model the atomic interaction by a 2D contact potential 
characterized by, 

Hrnt = ^Y.5'-^^^''^~''^) (6) 

i<j 

where g — ^/Sna/Xz is the dimensionless coupling, a is the 
3D scattering length and = y'h/Muj^. We assume 
the trap frequency in the z direction much larger than any of 
the energy scales involved, in such a way that only the lowest 
level is occupied, the dynamic of the system is frozen in the z 
axis and can be considered as two dimensional. 

The analytical solutions of the single-particle isotropic 
problem (t = 0) is given by the Fock-Darwin wave functions 

im0 / oTl 

where n (= 0, 1, 2, ...) is the Landau level, m (m > —n) is 
the single-particle angular momentum and is the Laguerre 
polynomial |32]. 

Within the lowest Landau level (LLL) approximation (n = 
0), we choose the set of Fock-Darwin functions given by, 
(/iOm(^, — ^ ^ I e~''^/^ r"' to represent operators and 
functions in the second quantized formalism. The single- 
particle eigenenergies are 

Eom = h{uj± — B*)m + hjj±. (8) 

To formulate the many-body problem, we consider the set 
of many -body Fock states |ni, ^2, ...) where rii refer to the 
occupation of the single-particle states (porm ■ From now on, 
we will omit the Landau level index n = Q. We truncate the 
single-particle state label m where rrimax and consequently 
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the dimension of the Hilbert space is fixed by the requirement 
of convergence of the resuhs. We perform exact diagonaliza- 
tion and for convenience, analyze the Hilbert space of many- 
body Fock states in subspaces with fixed total angular mo- 
mentum L. In general i is a non- well-defined parameter We 
want to stress that the convergency condition in much more 
demanding for the Hall response function than for the GS, as 
for the Hall response all the excited states are involved and 
this is especilly the case if some anisotropy is included. 



B. Optical lattice implementation 

We now propose an optical lattice implementation that pro- 
duces the rotating trapped system assumed for each site in the 
previous section. The spirit of the proposal follows the work 
of Popp et al. ll33ll . using an optical lattice potential to isolate 
a few atoms per site and simulate the gauge field through a 
fast rotation per site. Unlike the setup by Gemelke et al. 1.34.1 . 
our proposal relies on a optical superlattice with square geom- 
etry to create the rotation, not on the controlled interference of 
multiple beams with a time-averaged triangular lattice poten- 
tial. 

More precisely, we consider a trapping potential in the z 
direction stronger than any energy scale in our problem, so 
that our system can be regarded as purely bidimensional as 
was pointed out previously. We will assume a lattice poten- 
tial deep enough to supress tunneling between sites; in other 
words, each lattice site is to be regarded as an incoherent copy 
of the same experiment. The presence of the lattice is there- 
fore important not only for confining a few-particles system, 
but also for signal amplification. 

We will now explain how the trapping potential can be im- 
plemented by modulations of the laser intensities. The labo- 
ratory coordinate frame will be denoted with uppercase letters 
(X, Y), while the rotating frame will be denoted with lower- 
case (x, y), as in Sec. Ill Al These reference frames are related 
by a 2D-rotation matrix of angle fit. The harmonic expansion 
of each lattice site around its minimum yields the trap poten- 
tial described in the rotating frame in Eq. (|5)- This potential 
can be realized by three pairs of laser beams in a standing- 
wave configuration: two of them with the same wavelength A 
in the X and Y directions, respectively, and a third one with 
wavelength A' — \/2A in the tilted X+Y direction (i.e., along 
the line Y = X). The laser intensities associated to this con- 
figuration are (see the Appendix) 



I{X,t) =. Ix{t)sm\kX), 
I(Y,t) = lYit)sm^ikY), 
I{X,Y,t) = IxY it) sin' ik' (X + Y)), 



t = 



t = 7r/2n 



V(x,y)/Vo 



(9) 



where 



Ix{t) = V'o[l+4r(cos(2m)-sin(2m))], 
lyit) = yo[l -4r(cos(2m) +sin(2m))], 
Ixvit) = VqSt sm{2nt). 



• • I 




(10) 



FIG. 1 : Contour plot of the sum of potentials given by the intensities 
in Eq.9 with r = 0.1 at t = (left) and t = 7r/2r2 (right). The figure 
shows how the resulting potential is a square lattice of anisotropic 
harmonic oscillators which are rotating individually with frequency 
Q. d = n /k, k being the wavelength of the laser. 



The intensity modulations tune both the anisotropy r and the 
rotation frequency 17, while preserving the average trap fre- 
quency, Vb = 7\fa;^/(2fc^) . In particular, Eq.(19) in the Ap- 
pendix shows an explicit derivation of the trapping frequency 
uj± as a function of the parameters that define the lasers build- 
ing the optical lattice: 

ujI = —Re{a)k^El 

Figure 1 displays a contour plot of the periodic trapping 
potential for two different times. Deformation and rotation is 
shown explicitly. 

The state is prepared by adiabatically loading a BBC into 
the optical lattice with a fixed small anisotropy r and a pre- 
defined rotation Vl for each site. The system will be initially 
in the approximate ground state, where the chemical potential 
defines the inhomogeneous atom number density in the cen- 
tral region of the experiment. Once the system is in a ground 
state with possibly nonzero angular momenturn, the lattice site 
anisotropy r is adiabatically switched off 11331 [3511 . Following 
this procedure, the starting symmetric stationary state per site 
is prepared and its linear response can then be analyzed. 



m. LINEAR RESPONSE FUNCTION 

Once we have the possibility to simulate magnetic fields 
acting on charged particles in an effective way, in the next 
step we are going to characterize our interacting many-body 
states using its susceptibility under weak perturbations. The 
philosophy is to perturb the system using a small oscillating 
term that moves the system in one direction (say, the y direc- 
tion) generating mass current and to measure the response of 
some observable A that captures the torque experienced by the 
density distribution due to the presence of an effective mag- 
netic field perpendicular to the xy plane. Significant values of 
the observable A mean significant ability of the unperturbed 
system to manifest Hall effect. Inspired by a recently pub- 
lished experimental observation of the superfluid Hall effect 
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11251] . we made the appropriate selection of the observable A 
and the perturbation. 

We calculate the linear Hall response of the sequence of sta- 
tionary states generated at increasing values of f2. Our goal is 
to characterize the many-body states by quantifying thek Hall 
behaviour and analyze the role of the excitation spectrum in 
the response. Once the diagonalization is done for a fixed il, 
the eigenstates and eigenvalues of the Hamiltonian are known. 
Let us call them and E^, ly = 0,1, 2, .., — 1, respec- 
tively, riij being the dimension of the Hilbert space considered 
in the diagonalization. 

We choose for the periodic perturbation an extra term in 
the Hamiltonian H — Hq + Hpert given by Hpert — 
M euj\e^^cos{Lot)B where B = J2f and for A = 

AJ^ X]f (^y)i ' ^ being a dimensionless small parameter to 
ensure a perturbation treatment. Had we tried with the op- 
erator B ~ yi (equivalent to a constant force in the y 
direction), some nonzero Hall response would also have been 
obtained ||25I1 . However we would not expect any contribution 
from particle interactions, a necessary ingredient in the FQH 
effect, because this dipolar perturbation produces the center- 
of-mass displacement of the whole system. However, further 
analysis is necessary in this direction. 

The parameter ry is assumed to be small enough to en- 
sure adiabatic evolution of the system starting from to — > 
—oo when the perturbation is negligible, up to the stationary 
regime. From standard linear response theory at zero temper- 
ature Ii37i] we obtain 



A{A) = 2 e |x(w)| cos{ujt + 6) 



(11) 



where A (A) means the change of the expected value of A 
from the remote past when the perturbation was not active, 
to the moment when the measurement is performed. S is the 
phase of the complex (dimensionless) given by 



mw){i^\A\o) 

— + (jj + iri 



{0\AW){i^\B\0) 
Ely — Eq — uj — irj 



i<5(Lj) 



(12) 



The sum is extended to all excitations. However, only 
quadrupolar excitations have non-zero contributions due to 
the quadrupolar nature of the perturbation B. Being more ex- 
plicit, in the second quantized formalism the operators A and 
B take the form 



B 



+ \/(m + l)(m + 2)a]„a„i+2 — 2(m + 1) 



— E(v/m(m- l)a\^ar. 



^{m + l){m + 2)a]'„a,„+2 ) , 



(13) 



which can only change the ground-state angular momentum 
in two units. 



Within our lattice proposal, the small driven perturbation 
Hpprt can be implemented by a slight modification of the tun- 
able lattice intensities given by 



5Ix{t) = 2eVbcos(wt)[sin^(m) +sin(2m)/2], 
5lY{t) = 2eVocos(wt)[cos^(m) +sin(2m)/2], 
6IxY{t) = -2eVbcos(a;t)sin(21]i), (14) 



which produce the required modulated perturbation, as shown 
in the Appendix. This perturbation is maintained for a time 
T until the stationary-state regime is achieved. The lattice 
is then switched off abruptly and a time-of-flight (TOP) im- 
age of the system is taken. The fast expansion of the atoms 
maps the Fourier transform of their wavefunctions to position 
space. Since many copies of the experiment are performed 
at the same time, time-of-flight images provide the expecta- 
tion values of the momentum operators, (p^Py), which are 
fourier related to the observable, (xy), which we study nu- 
merically in the next section. It is important to remember 
that the TOP measurements take place in the laboratory refer- 
ence frame. Thus, (pxPy) will have to be reconstructed from 
the actual measurements through the relation {pxPy){t) ~ 
{PxPy){t) cos(2m)/2 + (P2 _ p2)(i) sin(2m)/2. This can 
be done by inverting the unitary rotation matrix or via filtering 
with a frequency Vl. 

Finally, we will address the experimental feasibility of our 
proposal. First we analyze the independent lattice sites ap- 
proximation. We can set an upper bound to the tunneling 
parameter as J < Vo5(Vo) where S is the overlap between 
the ground-state wavefunctions for neighbouring lattice sites. 
Let us assume a very deep lattice Vq > 'iQEn, where En = 
fi^ /2m}? is the recoil energy, which for '^^K and ^"^ Kb in a 
A = 800nm lattice is about Eji/h ^ 50kHz. For the ground- 
state wavefunction with L — A, chosen to match most of the 
results shown in the numerical simulations in section |IV] and 
considering isotropic lattice sites, we numerically estimate the 
tunneling parameter between neighbouring sites to be in the 
order of a few kHz. This value is just 10 times larger than 
when considering a ground-state with L = 0, which means 
that Vb can be kept constant at all times. Secondly, we must 
show that the quadrupolar excitations which can be probed 
with our Hall response perturbations are lower in energy than 
the lattice bandgap, so that the single-well experiment approx- 
imation is fulfilled. For a deep lattice of Vn = SOi?^, this 
bandgap can be estimated as Af,g ^ lOE'^j IBol . or At,g/h ^ 
300 kHz with the same choice of experimental parameters as 
before. We can see in Fig. |3] noted below, that in our proposed 
rotation regime i7 = 1.8 (1.8a;^/2 ^ 500 kHz) the first ex- 
cited states have an energy of about E„ — i?o ^ (0.1 — 0.3) 
or (0.1 - 0.3)a;±/2 - (30 - 80) kHz, so that additional bands 
will not be significantly populated. It is worth noting that, if 
O is further increased, and the regime for L = 12 is reached, 
the value of Vq must be reviewed. 
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1.74 1.76 1.78 1.80 1.82 1.84 1.86 1.1 

n 



1.90 1.92 1.94 



FIG. 2: Angular momentum of the ground state in units of h, as a 
function of Q in units of oj±/2 for several anisotropies. t — solid 
line; r = 0.8 x 10"^ dashed line; r = 5.0 x 10"^ dotted line; 
and T = 0.025 shown by large dots. The largest possible value for 
B* decreases as the anisotropy grows, being B* — 2AIQ = 1.897 
in our units for r = 0.025 




FIG. 3: Energy gap, Ei, — Eq versus fl for the first step from Q — 
1.76 to n = 1.92 with Las = 4 (see Fig.2). We show only the 
energies of the quadrupolar excitations, namely, those with L — 6 
(solid lines) and L = 2 (dashed lines). The horizontal dotted line 
at 0.3 selects a value for uj to visualize the crossings at four values 
of r2( see Fig.5(b) below). The energies are in units of hu}±/2, the 
frequencies in units of a;x/2, and the angular momenta in units of h. 



isotropic case (r = 0, full line), only some special values 
of the angular momentum are possible GS's, the so-called 
"magic numbers": L = 0, is the condensed state (with fi- 
delity = 1 , i.e., the overlap between the exact solution and 
the analytical expression), L = 4 is the Pfaffian state (with fi- 
delity = 0.95), L = 8 is the quasiparticle state (with fidelity 
= 0.98 ) and L = 12 the Laughlin state (with fidelity = 1 ) 
m. At the values of = 1.76, 1.92 and 1.96, transition 
jumps (steps) between different angular momenta take place, 
precisely where several eigenstates with different angular mo- 
mentum become degenerate. These are the only values of Q. 
where anisotropic configurations of the ground-state (e.g. , two 
voitices) are possible as linear combinations of the degenerate 
states. The analytic expression for the location of the first step 
is = cj^(l - gN/iSn)) 121, where N is the number of 
particles. Without loss of generality, we assume gN = 3, 
for which ilc = 1.76. If some anisotropy is considered 
(t = 0.8 X 10"^ or r = 5.0 x 10"^ in Fig.2), the steps are 
smeared out. The transition takes place in a continuous way, 
over ranges of of finite size within which anisotropic con- 
figurations are possible. If the anisotropy is large (t — 0.025 
in Fig.2) the step structure disappears and is replaced by an 
increasing monotonous function, as shown by the large dots. 
If instead we increase the number of particles maintaining a 
small anisotropy, the efect on the function {L)gs/^ is quali- 
tatively similar: The number of steps increases, a sequence of 
micro-plateaux appears [38J and for large values of 51, Lqs 
becomes a nearly continiuos increasing function similar to the 
one shown in Fig.2 for large anisotropy. 

Figure 3 displays the spectrum (i?^ — Eq) versus over the 
plateau Lqs = 4 in the isotropic case. Only the quadrupolar 
excitations to L = 6 and L = 2 are considered, since they 
are the significant ones for our choice of perturbation and ob- 
servable operators A and B (the operators and xy change 
L in ±2 units; see Eq.(13)). Once lj is fixed, every cross- 
ing of a constant horizontal line at uj with one of the lines of 
the spectrum is a candidate to be a peak of x{^) I ^ where, the 
resonant condition cancels the denominator in the second term 
of the right hand side of Eq(12). For example, if w = 0.3 , as 
shown in Fig. 3 there are four crossings between the horizon- 
tal line at 0.3 and the excitations with L ~ ^ (with negative 
slope) or L = 2 (positive slope), as shown in Fig. 5(b), or for 
UJ = 0.4 there are two crossings, as shown in Fig. 6. 



B. The Hall response 



IV. NUMERICAL ANALYSIS 

A. Ground and excited states 

From now on we will consider four particles, unless other- 
wise stated. Before showing the results of the Hall response 
function x{^)^ it is convenient to have in mind the results 
in Figs. 2 and 3. Figure 2 shows the angular momentum of 
the GS as a function of fl for different anisotropies. For the 



For the symmetric case. Fig. 4 shows a detailed analysis of 
the structure of close to the first step at ilc as uj grows 

from zero. The small peak (visible in Fig. 5(a) for a; = in 
a different scale) splits into two peaks which separate from 
each other. The important feature is that the height of the step 
at r^c, not related to any resonant condition, also depends on 
cj. The small peak moves to the light following the cross- 
ings between the positive-slope line of the lowest excitation 
with L — 2 (see Fig. 3) and the values of uj. In Fig. 5 four 
typical cases are shown for the whole range of starting at 
1.74 . Two different scales are considered: Figs. 5(a) and 5(d) 



6 




FIG. 4: Modulus of the Hall response x('^) (dimensionless) for dif- 
ferent values of cj. The tic-labels of the x-axis correspond to the low- 
est curve, the rest are shifted for clarity. The y-axis is also shifted. 
The arrows marc the first step at fixed flc ( = 1.76127), independent 
of ui. The frequencies are in units of cjx /2. 




1.76 1.80 1.84 1 



1.76 1.80 1.84 1.88 1.92 1.96 

n 



FIG. 5: Modulus of the Hall response x(.^) (dimensionless) for dif- 
ferent values of lj. From (a) to (d) uj — 0.0, 0.3, 0.5, 1.1 . The 
whole range of Q, values from 1.74, where the lowest Landau ap- 
proximation is still valid, to 1.98 close to the upper limit at 2 is con- 
sidered. The frequencies are in units of /2. 



up to 10 and Figs. 5(b) and 5(c) up to 100. The first plateau 
from n = 1.76 to 1.92 is fully dominated by the dynami- 
cal response peaks. In contrast, the dynamical response de- 
creases for larger L. For Lqs > 4, the peaks disappear be- 
yond u! — 0.4 . For uj — 1.1 , close to the driving frequency 
chosen in Ref. jzsll . the peaks completely disappear. 

Two comments related to Fig. 5 are in order: The physical 
expected behaviour of the response is given by finit values in 
the limit 77 ^ (see Eq.(12)). In our simple model, this 
is not the case as |x(w)| diverges in this limit at the resonant 
points. The inconvenience comes from the fact that we are not 



considering the width of the excitations that would prevent 
the divergence. The alternative procedure followed here is 
to consider a cutoff of such that we are close but out of 
resonance and the results are rj independent. Therefore we 
obtain information about the presence of a peak but not of 
the details close to its center, where non-linear response is 
expected due to absorption. At the end, we characterize the 
ground state. 

The other comment refers to the absence of peaks for val- 
ues of SI and cj where there is a crossing (Ei, ~ Eq = u (see 
Fig. 3)). The explanation of this possibility is, however, much 
technical. The absence of a peak means that < 0|A|i/ > 
and/or < 0|B|i^ > cancels, even though > is a quadrupo- 
lar excitation. Every eigenstate at fi has associated a specific 
set of non-zero single-particle occupations, say {/?;} , where 
Pi must not be an integer number. Here i labels the single- 
particle Fock-Darwin function . defined previously. If the 
set defining A\0 > (or B\0 >) has no coincidences with 
the set defining the excitation, then the numerator in Eq.(12) 
cancels. 

Furthermore, from the results shown in Fig. 5 we obtain the 
following useful information in the line of the characterization 
of the strongly correlated states. As SI increases, the entan- 
glement of the GS grows since the sp occupations Pi equal- 
ize, producing stationary states with internal structure far from 
condensation and mean-field description. The extreme case is 
the Laughlin state with a nearly equal distribution of occupa- 
tions, even for finite systems, showing large entanglement fl. 
This large entanglement has two consequencies: The Hall re- 
sponse is large (see Fig. 5(d)) and the GS is protected against 
absorptions at odds with the expected behaviour as there is a 
large amount of possibilities detected as crossings in the spec- 
trum (see Fig. 3). 

Figure 6 displays the case a; = 0.4 with a slight anisotropy 
(t — 0.8 X 10"'^). As expected, only the steps are mod- 
ified due to the lifting of the degeneracy produced by the 
anisotropy. However, the peaks remain unchanged. 

As previously noted in Sec. Ill, the guideline of our perfor- 
mance has close connection with the first experimental obser- 
vation of the Hall effect with atoms ||25l] . A brief explanation 
of the main ingredients of the experiment is as follows. The 
initial state is a large ( iV ~ 10^ ) strongly deformed cloud of 
bosons in a superfluid regime. By atom-laser coupling, they 
submit the system to an artificial magnetic field B perpen- 
dicular to the cloud. Next the system is perturbed along the x 
direction with a time-dependent modification of the trap po- 
tential given by 5U ^ x'^cos{ujt) equivalent to a force linear 
in X. Finally, they measure the time evolution of the second- 
order moment < xy > of the density. Their main result is the 
oscillation of < xy > /t if B or zero otherwise. This 
B -dependent correlation transport (in the x and y directions) 
is the Hall response. In their Fig. 4, the one with which we 
contrast our results, they show the amplitude of the oscilla- 
tions of < xy > as a function of B for a fixed u . The result 
is a monotonous increasing function that closely follows the 
superfluid hydrodynamic equations up to a point where B is 
strong enough to generate vortex states; from this point, ex- 
perimental points depart from the hydrodynamic prediction. 
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and Hall response has lower values than those predicted by 
the model. 

According to our results, for small systems and negligible 
deformation, classifying the states by their ability to manifest 
the Hall effect is equivalent to classifying the states by their 
angular momentum L, that is, the phase transitions he at the 
steps. Furthermore, following this similarity between Hall re- 
sponse and angular momentum we are naturally brought to 
an extrapolation: If we increase N (increasing the number of 
steps in Lgs/^ fss)]) and simultaneously add some deforma- 
tion, we expect for a Hall response a monotonous increasing 
function as the one shown by the large dots of Fig. 2, in close 
agreement with Ref . Ii25il as discussed in the above paragraph. 

To complete our comparison, we analyzed the role played 
by the vortex states for large values of _B . To this end, we 
moved to = 6 and analyzed this posibility. The sequence 
of Las for 6 atoms is: L = 0,6, 10, 12, 15, 20, 24, and 30. 
In the step produced by the change of L from 10 to 12 the cre- 
ation of a two-vortex state is possible. A very small anisotropy 
is enough to mix states of different angular momentum and fa- 
cilitate the precise numerical calculation of the value of Q, at 
the vortex state. Different Hall response, related to the gen- 
eral tendency, was expected for a vortex state due to its dif- 
ferent set of sp occupations {f3i] as compared with a state 
laying on a plateau. This difference is related to the more de- 
manding convergency condition in the numerical calculation 
for the vortex case. For a system laying on a plateau, a sin- 
gle occupation is dominant (the degree of condensation of the 
state is high) but in contrast two of the occupations play an 
important role in the vortex state [29]. Unexpectedly, our re- 
sult and the experiment goes in the opposite directions; Fig. 7 
shows a slight increase. It is difficult to follow the numerics 
to infer any difference when increasing N, or more impor- 
tantly, when increasing the deformation (up to a quasi-one- 
dimensional system which is the case in the experiment), two 
possible reasons of the discrepancy. 

Finally, a comment about the term "Hall response". We 
have followed the nomenclature used in the experimental 
work 1I25I1 which is an attempt to mimic the behavoiur of real 
charges under magnetic and electric fields. We believe that 
the observation of the torque of the density when the system is 
displaced in one direction, which has a significant value only 
in the case of non-zero B*, is related to the Hall effect and can 
be used to clasify the states. The absence of response would 
clasify the state outside the set of Qi/-type states. However, 
it must be noted that it is far from the simulation of the ap- 
propriate transport equation given by — (JxyEy where the 
conductivity a or its inverse, the resistivity, show the Hall ef- 
fect characterized by plateaux when analyzed as a function of 
the magnetic field. 



V. SUMMARY AND DISCUSSION 

In this work we have studied a system of interacting bosonic 
atoms under a strong effective magnetic field. Such sys- 
tems are known to exhibit a variety of phases, ranging from 
Bose-Einstein condensation to highly correlated states (Pfaf- 
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FIG. 6: Modulus of the Hall response xi'^) (dimensionless) versus 
Q, for uj — 0.4 with a slight anisotropy (r = 0.8 x 10"'^) has been 
considered. The frequencies are in units of io±/2. 
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FIG. 7: Expected value of the angular momentum (in units of h, 
upper part) and the modulus of the Hall response xi^) (dimension- 
less, lower part) as functions of for A'^ = 6. The upward peak at 
fl = 1.941 (signaled by an arrow) corresponds to a two-vortex state; 
a; = 1.1 is considered. The frequencies are in units of cjx/2. 



fian states, Laughling liquids, etc). While progress in experi- 
mental manipulation of ultracold atoms is advancing steadily, 
there is a huge need for new tools to probe and learn about 
the physics of strongly correlated states. Our work shows that 
the linear response can be a very useful method to get infor- 
mation about many-body atomic systems, both from the point 
of view of phase transitions, signalling the changes between 
different symmetries (the transition from L = to L = 
e.g., involves a broken parity symmetry at the single particle 
level |29]), and also for accurately analyzing the role played 
by the excitations. In what follows, we first summarize the 
whole experimental protocol, next briefly explain our results 
and finally, present our interpretations and comments. 

We propose an experimental setup using an optical lattice. 
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where several incoherent copies of a few-particle Hall state 
can be prepared. We have shown that the system can be pre- 
pared in a ground state with fixed angular momentum: A BEC 
is adiabatically loaded into an optical lattice with local defor- 
mation and rotation and once the system is in its ground state, 
the deformation is adiabaticaly switched off while keeping L 
constant. Next the system is perturbed by using laser intensity 
modulations at a particular frequency uj, and when the station- 
ary state is achieved, the lattice is switched off and a time-off- 
flight imaging of the cloud is performed. This measurement 
provides the density distribution required to build the density 
moment < xy > (see Eq.(ll)) and from it, the linear Hall 
response xi'^)- 

Our results can be summarized as follows: In the isotropic 
case, the phase transitions related to the Hall response take 
place at critical values of ft where changes in L occur in a 
stepwise well-known variation. Within a fixed phase, peaks 
at some specific values of ft and uj provide information about 
the excited states and characterize the dynamical response of 
the ground state within this phase. More importantly, we find 
a relationship between the Hall response and the correlation 
in the GS. As the correlation (or entanglement) increases, the 
Hall response grows; in contrast, the dynamical response ex- 
plicited by the resonant absorption peaks nearly disappears 
for perturbation frequencies comparable to the excitation en- 
ergies. For a vortex state of six atoms, we obtain a slight in- 
crease of the Hall response as compared with states without 
vortices at close il. 

Signatures of Hall response would be obtained if a stepwise 
structure of |x| is observed for large uj (see Fig. 5(d)), each of 
the plateaux related with a correlated ground state. One of the 
advantages of our implementation having isolated sites is that 
there is no influence of the lattice on the structure of the GS 
as is the case of other possibilities proposed to create FQH 
states, where a compromise must be achieved: The magnetic 
flux piercing each cell must be strong enough to produce the 
FQH state and at the same time, small enough to cancel the 
influence of the lattice avoiding the modification of the state 
ll39ll . The identification of the correlated states could be pos- 
sibly complemented by a local measurement of the two-body 
correlation function. Other ways to see the correlations be- 
tween the particles has been proposed in Ref. ll39[ l40ll using 
Bragg spectroscopy. 

The presence of an excited state is expected to produce a 
large response at a particular value of u! and consequently, 
energeticaly determined. In the case of vortex states, our re- 
sult, which is at odds with the behaviour shown in [25], is not 
conclusive. Extra analysis therefore is necessary for a larger 
number of atoms and/or vortices. However, having different 
Hall response as compared with the states without vortices, its 
presence can be detected experimentaly. Furthermore, accord- 
ing to our analysis, the stepwise variation of the Hall response, 
can be used to infer the angular momentum of the isotropic 
initial ground state. 

Two extensions of our analysis could, in principle, be ob- 
tained from experimental results: One is the inspection of the 
response close to the absorption peaks, inaccesible to our lin- 
ear response calculations. The other is the analysis of the evo- 



lution of the Hall response as N is simultaneously increased 
in each site, as a test of the expected extrapolated results. 
More importantly, the experimental difficulty that so far has 
prevented the observation of FQH states for large systems is 
expected to be solved in our setup, since lower rotation fre- 
quencies far from instability at are required to reach the 
FQH regime. Furthermore, from our optical lattice implemen- 
tation, we obtain an analytical expression for the trap potential 
frequency uj± as a function of the laser beams and the atoms 
involved in the experiment. 

The ideas in this manuscript are intimately related to the 
work by LeBlanc et al I25il . However, we must point out 
that the experimental results reported in that paper were ob- 
tained with strongly deformed (with a nearly cigar shape) 
Bose-Einstein condensates in the hydrodynamic regimes, with 
a large number of particles and a purely mean-field treatment. 
On the contrary, the states that we have studied in this work 
involve a small number of atoms in a symmetric trap, with a 
strong dominance of the interaction for large rotation speeds. 
Moreover, we have resolved the transition from vortex-free 
regime to the regime with vortices, going beyond the hydro- 
dynamical analysis present in that work. 

We are indebted with the referee and with M. Rizzi 
for all their suggestions and comments. H.P. and N.B. 
are greatly benefited from discussions with Maciej Lewen- 
stein. This work is partially sopported by the Spanish MFC 
through the FPU grant No.AP 2009-1761, the EU through the 
PROMISCE project, the Spanish MINCIN FIS2010-16185, 
the ConsoUder CPAN project CSD2007-00042 and the Gen- 
eralitat de Catalunya Program under contracts 2009SGR502 
and 2009SGR21. 



VI. APPENDIX 

Our goal in this Appendix is to obtain the relationship be- 
tween the Hamiltonian used in Section IIA for a single site 
(in the rotating frame of reference) and the appropriate laser 
beam implementation in the laboratory frame to reproduce it. 
To start, we rewrite Eq.(2) as 

Hsp — + Vtrap — ^L, (15) 

where Vtrap is given by Eq.(5) and L is the angular momen- 
tum operator in the z direction. From this expression, it is 
clear that the only term that must be translated from the ro- 

2 

tating to the laboratory frame is Vtrap since is invariant 
[41,1 . To this end, we analyze Eq.(5) in two terms: 

Vtrap{x,y) = Y(w^a;2 + w^y2^ 

= ^MujKx^ +y^) + 2tMujI{x^ - y^) 
= V1 + V2 (16) 

or following the convention given in the text about upper and 
lowercase letters, in the laboratory frame the trap potential 
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reads 

Vi{X,Y,t) ^Vok'iX^ + Y^), 

V2{X,Y,t) = 2tMuj\ X 

[X"^ cos(26') - cos(26') + 2XY sin(26l)], (17) 
where we have used the rotation matrix 

( x \ _ ( cos{0) sm{e) \ f X\ 
[y ) " l^-sin(0) cos(0) ) [y ) 

with e = nt and Vq = Mu}l/{2P) . 

We now consider on the one hand, that the intensity of 
a standing wave in the X , Y and X + Y directions are 
given by / = 8E^sm^{kX), I = 8E^sm^{kY) and / = 
8E^sm'^{k'{X + Y)) respectively (being k' = fc/y2), where 
Eo is the amplitud of the electric and k ( fc' ) is its wave vec- 
tor And on the other hand, that the coupling of the laser with 
the atomic induced dipolar moment d is d-E — 2Re{a)I{r) 
where / is the total intensity ifToll . Next, considering these two 
results, we can express Eqs.(A3) in terms of the laser-atom 
coupling in a periodic configuration. Previously, we want to 
stress that from Eqs.(A3) it can be inferred that considering 
only the first term in the expansion of sin^(fcX), two standing 
waves on the X and Y directions are sufficient to generate a 
symmetric time-independent trap potential, however, a third 
laser in the X + Y direction is necessary to deform and rotate 
it. We can rewrite the Eqs.(A3) as 

Vi{X,Y,t) - Vn{sin^{kX) + sin^{kY)) 

V2{X,Y,t) - 4TFo[sin^(fc^)(cos(26i) - sin(20)) 
- sin2(fcy)(sin(26l) + cos(26')) 
+2sin^{k'{X + Y))sm{2e)]. (18) 



We assumed that the atomic polarizability is aij = aSij 
and made the rotating wave approximation. Identifying terms, 
finally we obtain the main result 

32 

c^i = —Re{a)k'El (19) 

or in other words, we have expUcitly obtained the relation be- 
tween the trapping frequency and the experimental parame- 
ters of our configuration (the atomic polarizability a, the wave 
vector k, and the intensity of the lasers E^). Finally, the term 
of the Hamiltonian that generates a periodic Vtrap is given by 

Vtrap = Ix{t)sin^{kX) + lY{t)s\n^{kY) 

+IxY{t)sm^{k'{X + Y)), (20) 

where 

Ix{t) = yo[l + 4r(cos(2m)-sin(2m))], 

Wit) = Vo[l-4r(cos(20i) + sin(2m))], 

IxY{t) = Fo8Tsin(2m). (21) 
Similarly, the expression 

Hj,ert{y) = eMiol cos(a;^)y^ , (22) 

transformed to the laboratory frame reads 

HpertiX,Y,t) = eMw^cos(wt) X 

{X^ sm^{9) + Y^ cos^{9) - 2XY sm{e) cos{e)) 

~ 2eVo cos{ujt)[{sm'^{e) + i sin(26l)) sm^{kX) 

+(cos {0) + -sm{20))sm^{kY) 
-sin(26l)sin2(fc'(X + y))]. (23) 

In all the expressions, we considered only the first term in the 
expansion of sin^(fcX). 
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